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The interaction between surge and rotating stall in an axial flow compressor was investigated from the viewpoint 
of an unsteady inner flow structure. The aim of this study was to identify the key factor that determines the 
switching phenomenon of a surge cycle. The main feature of the tested compressor is a shock tube connected in 
series to the compressor outlet through a diaphragm, slits, and a concentric duplex pipe: this system allows surge 
and rotating stall to be generated by connecting the shock tube with the compressor, or enables the compression 
plane wave injection. The unsteady characteristics and the internal flow velocity fluctuations were measured in 
detail, and the stall cell structure was averaged and visualized along the movement of the operation point under a 
coexisting state of surge. A coefficient of the cell scale fluctuation was calculated using the result of the averaging, 
and it confirmed that the processes of inner flow structure change differed from each other according to the next 
cycle of the surge. The result suggests that the key factor that determines the next cycle is the transformation of 


the internal flow structure, particularly between the stall cell and the entire circumferential stall, in both the reco- 


vering and stalling processes. 
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Introduction 


The internal flow of turbomachinery is intrinsically 
unsteady, and the flow structure consists of complex and 
entangled unsteady phenomena of fluid mechanics. In 
particular, while the compressor operates at the partial 
discharge area, grasping the feature of the structure is diffi- 
cult because of this complexity. The surge and rotating 
stall are well known as notable unsteady phenomena that 
are generated in axial flow compressor systems, and an 
array of reports have been published on these topics [1]. 
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Surge is considered as a one-dimensional self-induced 
oscillation, which occurs in the axial-direction of the 
compressor, and results from the negative resistance of 
the entire compression system. Consequently, surge ge- 
nerates large-amplitude fluctuations in pressure and mass 
flow at a low frequency. For this reason, surge is often 
called “global instability” [2,3]. In contrast, rotating stall 
is an instability in the circumferential direction. The 
low-energy region is originated in a local flow separation 
on the blade surface at low flow rates, and rotates as a 
stall cell in a circumferential direction at approximately 
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Nomenclature 
Ac stagger angle (°) 
A Blockage rotor blades clearance (mm) 
D diameter (mm) 
N compressor rotational speed (min ') 
RByos rotor blades position 
Fa revolution time of rotor blade (s) 
Tp relative position between stall cell and rotor 
U rotor tip speed (m/s) (= nD;N / 60) 
V number of stator vanes 
Voi circumferential velocity (m/s) 
Wp wall pressure-rise (Pa) 
Z number of rotor blades 
Ps static pressure-rise (Pa) 
S flow passage cross section area 
t time (s) 
u time-averaged axial velocity (m/s) 


half of the rotor speed. Because of these characteristics, 
the rotating stall is often called a “local instability” [4, 5]. 
Surge and rotating stall have been regarded as completely 
different phenomena, owing to the differences in their 
characteristics and formation. The oscillating phenomena 
depend on the ratio of the Helmholtz resonator frequency 
to the frequency of the rotating stall. This ratio is called 
Greitzer’s B parameter, and the experimental results in- 
dicate that the critical value of the B parameter, which 
classifies surge and rotating stall, is approximately 0.8 
for an axial flow compressor [6]. 

Tryfonidis et al. [7] and Garnier et al. [8] found 
through detailed experiments and numerical simulations 
that rotating stall plays an important role as a precursory 
phenomenon to surge. On another front, McCaughan et al. 
[9, 10] examined the bifurcation phenomena between 
surge and rotating stall, and determined that the boundary 
between surge and rotating stall is not only affected by 
the B parameter, but also depends on the steady-state 
compressor and its resistance characteristics. Following 
this, research on the interactions between coexisting 
surge and rotating stall progressed to a Greitzer-Moore 
model analysis [11, 12], demonstrating the correlation 
between these phenomena and the total-to-static pres- 
sure-rise characteristics of a compressor. However, the 
interaction of instabilities in the coexisting states of surge 
and rotating stall has not yet been investigated in detail, 
from the point of view of the inner flow structure and 
unsteady behavior. 

Our previous study [13] has shown that the size of the 
rotating stall cell is closely linked to the unsteady beha- 
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Greek letters 

o flow coefficient 

Paty divided section of flow coefficient 
p atmospheric density (kg/m°) 

y pressure-rise coefficient 
Subscripts 


first stage 


2 second stage 

c casing 

h hub 

n natural number 

r rotor tip 

Abbreviations 

BPF band pass filter 

CBRF coefficient of blockage ratio fluctuation 
LPF low pass filter 

DPLAT double-phase-locked averaging technique 
CDPLAT conditional DPLAT 


vior of the surge cycle. In addition, the presence of an 
irregular surge cycle, in which different cycles are ran- 
domly selected, has been reported, and it has been noted 
that the internal flow structure may be an important fac- 
tor in determining the next surge cycle [14]. 

In the present work, a two-stage axial flow compressor, 
which has a shock tube equipped with a capacity tank 
connected in series to the compressor outlet, was used to 
investigate the various types of phenomena under the 
coexisting states of surge and rotating stall. In particular, 
the switching phenomenon of the irregular cycle was 
investigated in relation to inner flow structure. 

In the experiments, the unsteady characteristics of 
surge and internal flow velocity in each stage were 
measured and visualized using the double-phase-locked 
averaging technique. Furthermore, with the result of the 
averaging, a coefficient of blockage ratio fluctuation was 
calculated. It confirmed that the processes of inner flow 
structure change differed from each other according to 
the next surge cycle. The results indicate that the key 
factor that determines the next cycle is the transformation 
of the internal flow structure, particularly between the 
stall cell and the entire circumferential stall, in both the 
recovering and stalling processes. 


Experimental Procedure 


Experimental apparatus 

The tested compressor used in the experiments is 
schematically shown in Fig.l, and its design features are 
listed in Table 1. The compressor is an expandable-stage 
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Fig. 1 Tested axial flow compressor and measuring system. 


axial-flow compressor having one to three stages with 
variable rotor blades and stator vanes. The number of 
rotor blades Z is 12, and the setting angle of the blade can 
vary from 25 to —25° in increments of 5°. On the other 
hand, the number of the stator vanes V is 15, and the 
mounting angle of the vanes can be set independently 
from 30 to —30°. The rotational speed N is 12,000 min ', 
and the rotary valve set on the exit duct adjusts the oper- 
ating point. In these experiments, the tested compressor 
is set to a two-stage axial compressor specification, with 
the rotor blades and stator vanes set to the design mount- 
ing angles. 

The main feature of the tested compressor is a shock 
tube with a capacity tank of 0.262 m’, which is connected 
in series to the compressor outlet through a diaphragm, 
eight slits, and a concentric duplex pipe, as depicted in 
Fig.l. When the compressor is connected to the shock 


Table 1 Features of tested compressor 


Tested axial compressor 


Rotational speed N 12000 min! 
Airfoil configuration NACA65 
Hub diameter Dı 80 mm 
Casing diameter De 131 mm 
Rotor blade 

Number of blades Z 12 

Rotor tip diameter D, 130 mm 
Chord 30.0 mm 
Clearance 0.5 mm (tip) 
Stagger angle (at tip) 64.06° 
Stagger angle (at hub) 30.50° 


Stator vane 


Number of vanes 4 15 

Chord 24.5 mm 
Clearance 0.5 mm (hub) 
Stagger angle (at tip) 9.99° 

Stagger angle (at hub) 16.84° 


tube by removing the diaphragm, the B parameter in- 
creases over 0.8, thereby generating the coexisting states 
of surge and rotating stall. By contrast, when the com- 
pressor is disconnected, surge is not observed because the 
B parameter falls below 0.8. The area shown with the 
black arrows in Fig.1 is the flow passage of the discharge 
of the pipe, and the domain shown with the grey lines is 
the inner side of the duplex pipe connected with the tank. 
The material used for the diaphragm is rubber film. The 
compressed air stored in the capacity tank advances to 
the inner side of the duplex pipe by breaking the diaph- 
ragm with a plunge. It is fed into the compressor cascade 
through the slits as a compression plane wave. 


Measuring methods 


The steady-state compressor performance curve was 
determined using the data for the pressure rise characte- 
ristics and axial flow velocities. The data were measured 
with a static manometer and a one-dimensional single 
hot-wire anemometer (DANTEC, 55R01) installed on the 
outlet duct of the compressor. For this measurement, the 
diaphragm was installed in the outlet duct, and then the 
compressor was disconnected from the capacity tank. 

The unsteady pressure fluctuation was measured using 
a precision pressure transducer (Kulite, XCQ-093-25A), 
whereas the axial velocity was measured using a 
one-dimensional single hot-wire anemometer installed on 
the same point as the steady-state measurement system. 
The inner flow structure of the compressor was investi- 
gated via the internal flow and wall pressure fluctuation 
data. The flow fluctuation data were obtained by the 
one-dimensional single hot-wire anemometers, each in- 
stalled on a position 5mm ahead of the rotor blades. On 
the other hand, the pressure fluctuation data were meas- 
ured using a precision pressure transducer (Entran, 
EPIH-343-100P) set on the same position as the first 
stage internal flow measurement location. In order to 
investigate the surge behavior, the diaphragm was opened 
and the compressor was connected with the capacity tank. 
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The measuring system previously described was also 
used for the compression plane wave injection experi- 
ments. The photo interrupter was installed at the induc- 
tion motor as the blade trigger for the averaging proce- 
dure. In order to investigate the unsteady behavior of the 
stall cell during the surge cycle, the double-phase-locked 
average technique (DPLAT) was adopted. The details 
will be described later. 


Experimental Results and Discussion 


Stall characteristics and surge cycles 


The steady-state compressor performances are shown 
in Fig.2 (a) as black circles. The flow coefficient ¢ and 
the static pressure-rise coefficient y in the figure are de- 
fined as follows: 

ġ=u]U , y=2p,/pU? (1) 

The design operating point of the tested compressor, 
which is ¢= 0.360, is shown in Fig.2 (a) as a white cir- 
cle. According to the reduction of flow from the design 
operating point, the compressor operating point runs into 
two types of stall conditions, as shown in Fig.2 (a). The 
first step is considered as “mild stall,” in which only one 
stall cell generated near the rotor tip region rotates at 
approximately 55% of the rotor tip velocity in each stage. 
Furthermore, throttling the mass flow rate induced the 
second step of the stall, called “deep stall”, where the 
stall cell expands toward the rotor hub region and over 
the entire range of the tip region. 

The tested compressor was unsteadily operating on the 
negative slope region of the performance curve when the 
capacity tank was connected, and the coexisting pheno- 
mena of surge and rotating stall were generated. The un- 
steady transient performance curves of the surge cycles 
are shown in Fig.2 (b), (c) and (d) as solid lines in des- 
cending order of mass flow rate. The flow coefficient 
corresponding to the valve position at the steady-state 
operating point is shown as broken lines in the figure. A 
“large cycle” was observed from the neighborhood of the 
peak point of the performance curve to the flow rate of 
mild stall region. The unsteady operating points oscil- 
lated back and forth between the stable and deep stall 
regions during the “large cycles.” In Fig.2 (b), the flow 
coefficient is set to the higher flow region (¢ = 0.215), 
and then the unsteady flow coefficient recovers to the 
higher flow coefficient. In contrast, if the valve position 
is set to the partial flow rate region (¢ = 0.175), then the 
unsteady flow coefficient moves to the peak point of the 
performance curve without sufficient flow recovery, as 
shown in Fig.2 (c). According to the valve position, dif- 
ferent shapes of the surge cycles are observed. However, 
these cycles are, by definition, the same “large cycle.” A 
“small cycle” was observed by further throttling the mass 
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flow rate to under ø = 0.155, as seen in Fig. 2 (d). The 
unsteady operating points pass between the mild and 
deep stall regions. In the case of a single stage specifica- 
tion, a “stagnation stall” occurred after the size reduction 
of the cycle. However, in these experiments with a 
two-stage compressor specification, only a limited “small 
cycle” was observed. 
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Fig. 2 Stall and surge cycles of tested compressor. 
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Fig. 3 Surge cycle generated by external disturbance and 


hysteresis of compressor characteristics. 


On the other hand, a cycle caused by an external dis- 
turbance and the existence of hysteresis in the stable op- 
erating area is called a “top cycle,” where the unsteady 
operating points wander up and down between the stable 
operating region and the mild stall region. At the stable 
operating point ¢ = 0.245, a “top cycle” was observed 
caused by the compression plane wave injection, as seen 
in Fig.3 (a). On another front, when a surge cycle occurs 
on the negative slope of the performance curve, and even 
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if the valve was opened and it recovered to the stable 
operation point, the “top cycle” caused by the effect of 
hysteresis was also observed, as seen in Fig.3 (b). As 
stated above, a “top cycle” differs from the cycle occur- 
ring during the partial flow operation because it is ob- 
served in the steady state operation point. 

In addition, two types of irregular cycles were de- 
tected, in which two different surge cycles are selected 
irregularly. One is the “bottom irregular cycle” which 
occurs at the operation point of ¢ = 0.155, where the 
“large cycle” and the “small cycle” are randomly 
switched, as in Fig.4 (a). The other is the “top irregular 
cycle” at ø = 0.235, where the “top cycle” and the “large 
cycle” are replaced by each other in a chaotic pattern. 
The present study narrows the focus of the investigation 
on these two irregular cycles. 
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Fig. 4 Irregular surge cycles. 
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Conditional double-phase-locked averaging technique 


The internal flow field under the coexisting state of 
surge and rotating stall is intricate in construction. This 
complexity is attributed to the overlapping of the rotor 
blade rotational frequency and the fluctuating rotational 
frequency of the stall cell, which is associated with shift- 
ing the unsteady operation point in the blade cascade. 
Therefore, it is difficult to grasp the details of the stall 
cell features by using the double-phase-locked averaging 
technique [15], which is used to visualize a flow structure 
in the blade cascade under the assumption that the rota- 
tional speed of the stall cell is constant. 

In this work, in order to investigate the key factor that 
determines the next surge cycle of the irregular cycles in 
the coexisting state of surge and rotating stall, the inner 
flow structure was visualized by a conditional double- 
phase-locked averaging technique (CDPLAT). The rota- 
tional speed of the stall cell in the irregular cycles has a 
bandwidth of 90 to 130 Hz. Therefore, the band-pass 
filtered pressure fluctuation data were used as the vortex 
trigger, with the bandwidth set to the frequency range of 
90 to 130 Hz. 

For the preliminary preparation of the CDPLAT, the 
cycle change of each irregular cycle is classified into four 
cases. The bottom irregular cycle is classified into large 
to large (L-L), large to small (S-L), small to large (L-S), 
and small to small (S-S). The top irregular cycle is classi- 
fied into top to top (T-T), top to large (T-L), large to top 
(L-T), and large to large (L-L). 
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Fig.5 Procedure of conditional double-phase-locked averaging technique. 
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As an example, the time history of the experimental 
data of the “bottom irregular cycle” is shown in Fig.5. 
The figure shows the low-pass filtered pressure-rise coef- 
ficient y, the low-pass filtered flow coefficient ø, the 
circumferential velocity ahead of each rotor blades Vai ; 
and V.j,2, the band-pass wall pressure Wp (BPF) and the 
signal of the photo interrupter which locates the position 
of rotor blade. The workflow of the CDPLAT is as fol- 
lows. First, a reference point was defined as the local 
maximum point of Wp (BPF). Second, the data of each 
Vi, were delimited, and the length of the data from one 
reference point to next point was regarded as the length 
of one stall cell rotation. Third, the data were labeled 
based on the classification of the surge cycle change and 
flow coefficient section yn. Next, the average data 
length of one stall cell rotation was calculated based on 
the classification. Then, the data over 8 % far from the 
average were eliminated. After this filtering process, the 
remaining data were labeled based on their relative posi- 
tion between the stall cell and rotor blade. A pitch of the 
rotor blades was divided into eight regular sections, 
which had values Tp = 0-7. Once categorized, the data 
were averaged by category and then rearranged in accor- 
dance to their relative position between the rotor blade 
and the stall cell. 


Standardization and structure of stall cell in two 
stages 


As a preliminary step toward its visualization under 
coexisting phenomena, the stall cell occurring at the mild 
stall region of the steady state was visualized. Fig.6 
shows the circumferential velocity distribution ahead of 
each rotor blade. The left contour of the figure is the dis- 
tribution in the first stage, while the right one corres- 
ponds to the second stage. The circumferential mounting 
locations of the hot-wire anemometer and pressure 
transducer are marked in the figure. The influence of the 
stall cell and the wake of the rotor blades can be clearly 
seen. The configuration of the stall cell on its shape is 
similar to the previous reports by Das et al. [16] and 
Poensgen et al. [17]. As shown in the frame of the figure, 
the stall cell in the second stage precedes by about 30° in 
the direction of propagation as compared to the stall cell 
in the first stage. At the position of 3 o’clock in the first 
stage, an increase in velocity on the suction side of the 
blade can be seen, which may be related to the flow se- 
paration on the suction surface. On the other hand, at the 
position of 11 o’clock in the first stage, an increase on the 
pressure surface can also be seen. This increase may have 
relation to the flow that comes around the pressure sur- 
face from the suction side due to the separation point of 
the trading area of the cell. In comparison with the first 
stage, the overall circumferential velocity in the second 
stage 1s higher. 
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Fig.6 Visualization of the rotor inner flow structure using the 
double-phase-locked averaging technique. 


Investigation of stalling and recovering process by 
conditional double-phase-locked averaging technique 


Fig.7 shows the scope of application of the CDPLAT 
in these experiments. The switching phenomenon of the 
irregular cycle was observed in the recovering process of 
the bottom irregular cycle, as well as in the stalling 
process of the top irregular cycle. Therefore, in this study, 
the application range of the CDPLAT corresponding to 
the operation points where the cycle change had not oc- 
curred on the transient characteristics cycle of surge. The 
recovering process of the bottom irregular cycle was de- 
cided to be from the flow coefficient of ¢ = 0.157 (@j, 7) 
~0.252 (iy 9) divided into nine sections (Fig.7 (a)), while 
the stalling process of the top irregular cycle was to be 
from ø =0.232 (Qj, 1) ~0.144 (Pav 14) into fourteen sec- 
tions (Fig.7 (b)). 
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Fig. 7 Scope of application for the conditional double-phase- 
locked averaging technique. 


The distribution of the circumferential velocity ahead 
of each rotor blade is shown in Fig.8 for the recovering 
process of the bottom irregular cycle, and in Fig. 9 for 
the stalling process of the top irregular cycle. From top to 
bottom on the figure, the movement of the operating 
points are ordered, and the distributions are piled up for 
the cycles that have an identical last cycle. 
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Fig. 8 Distribution of the circumferential velocity ahead of 
each rotor blade in the bottom irregular cycle. 
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Fig. 9 Distribution of the circumferential velocity ahead of 
each rotor blade in the top irregular cycle. 


The decline of scale of the stall cell according to the 
recovery of the flow rate in Fig.8 was revealed, as well as 
the incline of scale following the stalling in Fig.9. In 
Fig.8, Qiiv2~ Qay sis the flow coefficient range discussed 
later, where the entire circumferential stall transforms 
into the rotating stall. Focusing on the outer frame region 
where the circumferential velocity is higher in the distri- 
bution, the decline in scale of the stall cell in the cycle 
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getting to the large cycle is faster than that getting to the 
small cycle. In particular, in the second stage framed with 
the solid line, a large decrease of scale was observed 
during Qj, 3to dj, 4. This decline preceded that of the first 
stage, Qi 4tO Qays. In Fig.9, in contrast, attention is paid 
to the enclosed dark color region, where the circumferen- 
tial velocity is lower. A growth in scale of the stall cell 
could be seen, and it is shown that the growth of the stall 
cell associated with the cycle getting to a large cycle 
preceded that of the top cycle. 


Coefficient of blockage ratio fluctuation 


In order to evaluate and compare the result of the 
CDPLAT quantitatively, a coefficient of blockage ratio 
fluctuation of the n-th is defined as follows: 


CBRF, -( j7 din ) -Vair ( mas [fas)r00 (2) 


U 


Using the CBRF calculation, the relative changes in 
scale of the stall cell were obtained. Fig.10 shows the 
result of the evaluation. 


From topcycle (b) 


From large cycle (a) 


From small cycle (c) From largecycle (d) 
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Fig. 10 Coefficient of blockage ratio fluctuation. 


The figures were displayed in a batch for the last cycle. 
The black circle frame shows the divided sections of 
flow coefficient where the bar charts were branched out. 
In the case of the bottom irregular cycle, the branching 
on the second stage preceded the first stage. On the other 
hand, in the top irregular cycle, a positive branching was 
shown at the first stage of each cycle, whereas a luke- 
warm branching can be seen at the second stage. In 
Fig.11, the areas enclosed in black circles were shown in 
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the characteristic curve as a yellow shading. From the 
figure, it was found that the branching on the bar chart of 
the CBRF calculation precedes the bifurcation on the 
unsteady characteristic curve. 

In Fig.12, in order to verify the tendency of the bran- 
ching from the viewpoint of raw data, the raw waveform 
data of the circumferential velocity ahead of each rotor 
blade were compared in relation to every cycle change. 
The shaded area falls under the category of both the en- 
closed area in Fig. 10 and the shaded area in Fig.11. 

In Fig.12 (al), comparing L-L cycle to L-S cycle in 
the first stage, in the recovering process of the bottom 
irregular cycle, the two waveforms tread a similar path in 
t/T,-y = 0~6, where the entire circumferential stall region 
transforms into the stall cell, and then reaches the mild 
stall region. While at the same time in the second stage, 


Fig.12 (a2), there is variability among the two waveforms. 


Furthermore, this trend of the pattern holds true when 
comparing the S-L and S-S cycles. It is reasonable to 
consider that this variability, observed during the trans- 
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formation of the stall region between the deep and mild 
stalls at the second stage, appear as the branching of the 
CDPLAT bar chart. The tendencies indicate the possibil- 
ity that the transformation process of the stall region 
ahead of the second stage rotor blade is the key factor 
that determines the next cycle during the recovering 
process of the irregular cycle. 
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Fig. 11 The region of operation points where the coefficient 
of blockage ratio fluctuation branches off. 
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Fig. 12 The raw waveform data of circumferential velocity ahead of each rotor blade. 
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On the other hand, in the stalling process of the top ir- 
regular cycle, the waveform has a different propensity to 
the recovering process. In Fig.12 (b1) and (d1), the two 
waveforms, which showed that the same fluctuation has a 
relationship to the rotation stall at first, gradually under- 
went a clear change. The cycle getting to the top cycle 
was continuously the same in waveform, whereas the 
cycle reaching the large cycle showed the form asso- 
ciated with a deep stall. This trend of change is applica- 
ble to those of the second stage, as shown in Fig.12 (b2) 
and (d2). However, the amplitude of fluctuation in the 
waveform of the cycle getting to a large cycle is larger 
than that at the first stage. The shaded area in Fig.12 (b1), 
(b2), (d1), and (d2) falls under the category of the region 
where the rotating stall transforms into the entire cir- 
cumferential stall. 

From the above, both for the recovering process of the 
bottom irregular cycle and the stalling process of the top 
irregular cycle, the results indicate that the key factor that 
determines the next cycle is the process of transformation 
between the rotating stall cell and the entire circumferen- 
tial stall. 


Conclusions 


The unsteady behavior of a two-stage axial compres- 
sor under a coexisting state of surge and rotating stall, 
and in particular the recovering process of the bottom 
irregular cycle as well as the stalling process of the top 
irregular cycle were experimentally investigated by de- 
tailed measurements of the unsteady characteristics and 
the circumferential flow velocity fluctuations ahead of 
each stage of the rotor blades. The internal flow structure 
during the irregular surge was visualized by using a con- 
ditional double-phase-locked averaging technique. Con- 
sequently, the increase and decrease in scale of the stall 
cell associated with the movement of the operation point 
could be seen. In addition, the calculation of a coefficient 
of blockage ratio fluctuation indicated that the key factor 
that determine the next surge cycle in a coexisting state 
of surge is the process of transformation between the 
rotating stall and the entire circumferential stall. 
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